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In the cold dark matter paradigm, structures form hierarchically, implying that large structures 
contain smaller substructures. These subhalos will enhance signatures of dark matter annihilation 
such as gamma rays. In the literature, typical estimates of this boost factor assume a concentration- 
mass relation for field halos, to calculate the luminosity of subhalos. However, since subhalos 
accreted in the gravitational potential of their host loose mass through tidal stripping and dynamical 
friction, they have a quite characteristic density profile, different from that of the field halos of 
the same mass. In this work, we quantify the effect of tidal stripping on the boost factor, by 
developing a semi-analytic model that combines mass-accretion history of both the host and subhalos 
as well as subhalo accretion rates. We find that when subhalo luminosities are treated consistently, 
the boost factor increases by a factor 2-5, compared to the typical calculation assuming a field- 
halo concentration. This holds for host halos ranging from sub-galaxy to cluster masses and is 
independent of the subhalo mass function or specific concentration-mass relation. The results are 
particularly relevant for indirect dark matter searches in the extragalactic gamma-ray sky. 


I. INTRODUCTION 

If dark matter is made of weakly interacting massive 
particles, their self-annihilation may produce high-energy 
gamma rays jlj]. Indirect searches for dark matter annihi¬ 
lation with gamma-ray telescopes are one of the promis¬ 
ing probes of non-gravitational interactions of dark mat¬ 
ter. In hierarchical structure formation, small structures 
form first and they merge into larger dark matter halos. 
Numerical simulations show that the distribution of dark 
matter particles in the halo is clumpy, with a substan¬ 
tial fraction being locked into substructures Hi- Since 
the self-annihilation rate depends on dark matter den¬ 
sity squared, presence of these subhalos will boost the 
gamma-ray signal. 

There are two well-adopted methods to estimate the 
boost factor Q. The first is to phenomenologically ex¬ 
trapolate subhalo properties, i.e., power-law scaling rela¬ 
tions between subhalos with a mass above a given thresh¬ 
old and their total luminosity, down to scales of the small¬ 
est subhalos (typically assumed to be on the order of 
Earth mass, although very sensitive to the exact par¬ 
ticle physics model U)> e.g., This approach yields 
very large boosts, on the order of 10 2 (10 3 ) for galaxy 
(cluster) halos, but there is no guarantee that this phe¬ 
nomenological extrapolation over many orders of magni¬ 
tude is still valid. In fact, this method is similar to a 
power-law extrapolation of the so-called concentration- 
mass relation. The second one relies on a concentration- 
mass relation that flattens toward lower masses. This 
behavior is favored analytically as well as from dedi¬ 
cated simulations 0 - 0 - Studies following this approach 
(e.g., I, 0, El) typically conclude that the boost fac¬ 
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tors are much more modest, about an order-of-magnitude 
below the phenomenological extrapolations Q 

The latter method is believed to yield more realis¬ 
tic values for the boost factor due to subhalos (B s h). 
The boost is typically calculated as an integral of 
('dN/dm)L s h(m ) over the subhalo mass m , with dN/dm 
the subhalo mass function as found in simulations and ex¬ 
trapolated down to the the minimal subhalo mass, and 
-bsh(ui) the subhalo luminosity, which is a function of 
the concentration. However, as mentioned by Ref. El, 
so far this method has not been used fully consistently, 
since the concentration-mass relation that goes into the 
calculation of L s h(m) is that of field halos, which is not 
directly applicable to the subhalos. In the gravitational 
potential of its host halo, a subhalo is subject to mass loss 
by a tidal force, which tends to strip particles from outer 
regions of the subhalo EHT5J. This effect will reduce 
the subhalo mass substantially, but keeps the annihila¬ 
tion rate almost unchanged, because the latter happens 
in the dense central regions dominantly. Consequently, 
subhalos are expected to be denser and more luminous 
than halos of equal mass in the field, and thus, the boost 
should be larger. 

In this paper, by developing semi-analytic models, we 
investigate the effect of tidal stripping of subhalos and 
show that using the field halo concentration indeed re¬ 
sults in a significant underestimation of the subhalo lu¬ 
minosity, and hence the annihilation boost factor. There¬ 
fore, we argue that this effect is extremely important in 
this context, and a consistent treatment of subhalo con¬ 
centrations should always be adopted. 

We note that there are alternative estimates of the 
boost factor that do not depend on the concentration- 
mass relation directly. Reference [l6i] applies an analytic 


1 See Appendix[C]for a calculation of the overall boost factor using 
the two different methods. 




2 


model for the probability distribution function of the halo 
density field including substructure. Reference sa uses 
a technique based on the stable clustering hypothesis and 
includes the effects of tidal disruption. Both of these are 
then matched to numerical simulations above the reso¬ 
lution scale. Finally, reference |l8] uses the nonlinear 
power spectrum directly to calculate the so-called flux 
multiplier, which encapsulates the boost, and thereby the 
extragalactic dark matter annihilation flux. 

We adopt cosmological parameters from 5-year WMAP 
results [l9| . Capital M refers to the host halo mass and 
lower-case m to subhalo mass. Quantities at redshift z = 
0 are denoted by subscript 0. Virial radius, r v ir, is defined 
as the radius within which the average density of a halo 
is A c (z)p c (z), where A c is given by Ref. [2C| and p c {z) 
is the critical density at redshift 2 . The virial mass is 
defined correspondingly. 


II. DENSITY PROFILE AND GAMMA-RAY 
LUMINOSITY 

The total gamma-ray luminosity of a dark matter halo 
of mass M is (e.g., [2l| ') 

L(M) = [1 + B sh (M)]L host (M), ( 1 ) 

B s h(M) = — —3_— / dm—L sh (m)[l + B ssh (m)], 

Thost(M) J dm 

( 2 ) 

with m the subhalo mass, B s h(M) the boost factor due 
to subhalos, Lh ost (M) and L s h(m) the luminosities of the 
smooth component of the host halo and subhalos, respec¬ 
tively (both often parameterized bv the Navarro-Frenk- 
White (NFW) or Einasto profile (22, El). According to 
the state-of-the-art numerical simulations, the subhalo 
mass function (i.e., number of subhalos per unit mass in¬ 
terval) behaves as a power-law dN/dm oc m ~ a , where 
a = 1.9-2, down to resolution scales [HI, 0, 25j. The 
boost factor due to “sub-substructure” B ss h(m) is either 
parametrized the same way as B s h or often neglected. 

Assuming that the density profile of the subhalos is 
characterized by the NFW function up to tidal radius r t 
(beyond which all dark matter particles are completely 
stripped), the subhalo luminosity is given by L s h(m) oc 
p 2 r 3 [1 — l/(l+c t ) 3 ], where p s and r s are the characteristic 
density and scale radius of the NFW profile, and Ct = 
rt/r s . In the literature where the effect of tidal stripping 
is ignored, one adopts the virial radius r v j r and virial 
concentration parameter c v i r = r v i r /r s instead of r t and 
Ct, respectively. 

III. ORDER-OF-MAGNITUDE ESTIMATE 

We start with an order-of-magnitude estimate. Rather 
than using physics-driven models, we rely on phenomeno¬ 
logical relations found in numerical simulations, where 


the effect of tidal stripping is automatically taken into 
account. Let us define Umax and r max as the max¬ 
imum circular velocity and radius where the velocity 
reaches Umax- For the NFW profile, these quantities 
are related to p s and r s through r s = r max /2.163 and 
p s = (4.625/47rG)(V max /r s ) 2 . 

For field halos, we assume the concentration-mass re¬ 
lation from Ref. J26| that matches well the simulation 
results of Ref. [T5|] down to the resolution limit. For 
a field halo of mass Wfh = 1O 5 M 0 , we find c v i r ~ 63. 
All other relevant quantities (r v n, r s , and p s ) then fol¬ 
low from mfh = 47rA Ci0 p Cj0 r 3 ir /3, r s = r vir /c v i r , and p s = 
TOfh/[47rrf/(cvir)], where /(c) = ln(l+c)—c/(l+c). From 
these, we find r ma x,fh ~ 40 pc and Vmax.fh « 1.2 km s . 

For subhalos, numerical simulations fl5[ found the 
following relation down to 10 5 M 0 : m s h ~ 3.37 x 
lO'(Vmax,sh/lOkms -1 ) 3 - 49 M 0 , from which we obtain 

Vmax.sh ~ 1-9 km s _1 for m s h = 1O 5 M 0 . The same sim¬ 
ulations found the relation between Vmax and r max for 
subhalos and those for field halos: (r m ax,sh/r m ax,fh) ~ 
0.62(V maXiSh /V max , fh ) 1 - 49 EH. Combining this with the 
results above for field halos and subhalos of equal mass 
(mfh = m sh = 10 5 M 0 ), we have r maXiSh « 50 pc. 

The ratio of the gamma-ray luminosity of 
the subhalo and field halo of mass 10 s M 0 is 
then L s h/Lfh ~ (Ps,sh/Ps,fh) (^Ysh/^s.fh/ 

(Vmax,sh/Vmax,fh) (Unax,fh/?"max,sh) ~ 5. We find 

that the luminosity ratio is weakly dependent on the 
mass. For example, A s h/Afh ~ 4 for m = 1O 9 M 0 . 
This result also holds for an Einasto profile, although 
there are some subtleties involved. For a more detailed 
discussion see Appendix ITU 


IV. SEMI-ANALYTIC MODEL 

Stripped subhalos tend to be denser than field halos of 
equal mass, and consequently more luminous. Below we 
quantify this difference in luminosity, which essentially 
depends on three parameters: r s , p s and Ct, all of which 
depend on the halo formation time, the infall mass and 
subhalo’s history in the host. 

We assume a truncated NFW function for tidally- 
stripped subhalos: p(r) = p s r 3 /[r(r + r s ) 2 ] for r < rt 
and 0 otherwise, in agreement with what is found in 
simulations 15(. Concerning the scale density and ra¬ 
dius, Refs. [27fl find from N-body simulations that the 
change in Vmax and r max , and consequently in p s and 
r s , only depend on the total mass lost by the sub¬ 
halo, following Vmax, 0 /Vmax, a = 2 0 - 4 a; cl - 3 /(l + x)°- 4 and 
r m ax, o/^max, a = 2 _0 - 3 a; 0 - 4 /(l+a;) -0 ' 3 , where x = m a /m 0 
and the subscript a represents epoch of accretion. 

Based on the extended Press-Schechter (EPS) formal¬ 
ism [ 2 H , Ref. dH provides an analytic model for the dis¬ 
tribution of infall times of subhalo progenitors into their 
host: d 2 A r /dlnm a /dln (1 + z a ) as a function of redshift 
z and host mass M(z). For the mass-accretion history of 
the host, we adopt the analytic EPS model from Ref. [301 ]. 
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This model provides the mean evolution of a halo that 
ends up with mass M 0 a,t z = 0. Therefore, we can pa¬ 
rameterize its mass at earlier times through M(z\Mq). 
Last, to take into account the effect of tidal stripping in 
the host we apply the semi-analytic model of Ref. [3l| . 
It provides an orbit-averaged mass-loss rate for subhalos, 
m (z\z a , m a , Mo). In this study, we assume their model, 
in which the mass-loss rate m/m is based only on the 
mass ratio m/M and the dynamical time scale, is valid 
for all mass-scales down to the smallest halos. This is an 
assumption that needs further testing and is the subject 
of future work. 

We start from a given set of two parameters that char¬ 
acterize subhalos, mo and c t . We solve the differen¬ 
tial equation for m(z ) backward in time in the gravi¬ 
tational potential of a host of mass M(z). Using the 
above-mentioned relations for p s and r s in terms of mo 
and m(z), we can compute the change in the tidal ra¬ 
dius and thus ct(z). For each step, we also compute the 
concentration-mass relation for the virialized field halo 
c v ir(m, z) Q, and once the subhalo Ct(z)-m(z) relation is 
found consistent with c v i r , we assume that the subhalo 
accreted at that particular redshift z a just after its viri- 
alization, and m{z a ) = m v i T (z a ) = m a . At this accretion 
redshift z a , the virial radius r v ; r of the subhalo is obtained 
by solving m v i r (z a ) = 47rA c (z)p c (z)r^ ir /3. The charac¬ 
teristic density and scale radius at accretion then follow 
from r Sia = r vir /c v i r and p Sja = m vir /[47rr 3 a /(c v i r )]- If 
the virialization happened earlier than z a , we would ob¬ 
tain a higher characteristic density p s ; therefore, our as¬ 
sumption is conservative. 

Finally, using these relations for m a and z a as func¬ 
tions of mo and Ct, and using the infall distribution, we 
compute a joint distribution function of mo and Ct: 


V (m 0 ,c t |M 0 ) oc 


d 2 N 

dmodct 


d 2 N 


dlnm a dln (1 + z a ) 
d (In m a , In (1 + z a )) 


d (mo, c t ) 


( 3 ) 


V. RESULTS 

We obtain our evolved subhalo mass function by inte¬ 
grating Eq. m over Ct, 

dN_ = f°° d 2 N 
dm 0 J Cmin dm 0 dc t Ct ' 

We take c m i n = 1 as the absolute minimum [32|. We 
checked that our results are insensitive to the exact choice 
as most halos have higher Ct- The maximum possi¬ 
ble value corresponds to the concentration of halos that 
formed and are accreted today. 


TABLE I. Properties of the evolved subhalo mass functions, 
dN/dmo oc mg a resulting from our analysis. Columns show 
the host halo mass, the mass fraction in subhalos assuming 
m m i n = 1(T 6 Mq, and the slope of the mass function. 


Afhost/MQ / su b 

a 

10° 

0.06 

1.93 

to 9 

0.08 

1.94 

10 12 

0.13 

1.94 

10 15 

0.23 

1.92 


Table Q] shows the characteristics of the subhalo mass 
function for host halos of different mass. The second 
column contains the total mass fraction in subhalos f su b' 

dN 

mo- dmo = /subAfhost^ (5) 

dm-o 

for which we adopted m m i n = 10 -6 Mq and m max = 
O.lMhost- The third column shows the slope of the mass 
function, a = — dln(d-/V/dmo)/dlnmo. We find good 
agreement with numerical simulations, only / su b being 
slightly lower (e.g., 0 ), 

We then compute the mean luminosity of a subhalo 
with mass mo = m, 

/ Cm ax 

Lsh(m, Ct)V (m, Ct) dct, (6) 

where T s h(m, c t ) oc p 2 r 3 [1 — 1/(1 + c t ) 3 ]. Figure Q] shows 
the luminosity-weighted mass function for subhalos in a 
Milky-Way-sized halo. Although the dependence is weak, 
smaller subhalos contribute more to the total subhalo 
luminosity. The upturn at the high-mass end is a result 
of the fact that the most massive subhalos can only be 
accreted at late times. Consequently, the evolved subhalo 
mass function looks more like the unevolved one, which 
has a harder slope. 



A. Boost ratio 

It is interesting to compare luminosities of subhalos ob¬ 
tained above with those of field halos of equal mass. We 
assume field halos to be virialized at z = 0 with r v i r given 
by m = 47rA c ,o/9 Cj oU) ir /3. The characteristic density and 
scale radius are again obtained with r Sj fh = rvi r /c v ir(m, 0) 
and p Si fh = m v ir/[47rr 3 /(c v ir(m, 0))], where the concen¬ 
tration mass relation of Ref. [9jj is assumed. Then the 
field halo luminosity is, Lfh(m) oc p 2 r 3 [l — 1/(1 + c v j r ) 3 ]. 
The dotted curve in Fig. [T] shows the luminosities Lfh(m) 
weighted by the same mass function as in the case of 
L s h(m). As anticipated above, the field halos are less 
bright than the subhalos of the same mass by a factor of 
a few, almost independent of mass 


2 See Appendix [3 for an alternative comparison using the Vinax 
rmax relation. 
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FIG. 1. Luminosity-weighted mass function for subhalos 
(solid), field halos whose concentration is set at 2 = 0 (dotted) 
and field halos which have the same infall times as the solid 
line (dashed) in a 10 12 M@ host. It shows the contribution of 
different mass subhalos to the overall subhalo luminosity. We 
always use the subhalo mass function from Table [I] 

It should be noted that subhalo concentrations depend 
on formation time. Halos that formed earlier are more 
concentrated since they formed in a denser background, 
an effect that has been taken into account in past studies 
(e.g., Refs. DlEHl). Since we set the concentration 
of the stripped halos at z a we also include the dashed 
line for a fully fair comparison. It shows the luminosity 
of halos that follow the same infall distribution as the 
solid line, and thus have the same natal concentration as 
this is set at the time of accretion, but are not tidally 
stripped. As can be seen, the tidal stripping still yields 
an increase by a factor of ~2 at any subhalo mass. The 
decrease in the difference in luminosity at lower masses 
is due to the smallest halos being accreted earlier, thus 
their concentrations at accretion differ most compared to 
that at z = 0. 

Since the boost depends critically on the subhalo mass 
function, in addition to our fully self-consistent model 
with the mass-function from Table Q] we also investigate 
dependence on several models for the mass function. We 
adopt four models, taking spectral indices of a = 1.9 
and 2, and smallest subhalo masses of m m ; n = 1O _6 M 0 
and 10 4 Mq. We compare the subhalo boost B s h(M), 
calculated with Eq. using subhalo luminosities of 
stripped halos, to the boost calculated without account¬ 
ing for tidal effects (using the virialized field halo mod¬ 
els) . Figure [2] shows the ratio of boosts as a function of 
host halo mass for these models. Taking tidal effects into 
account will enhance the boost by up to a factor of 5 com¬ 
pared to the simple field halo approach, consistently for 
host halo masses between 10 6 -10 15 Mq. This is largely 
independent of models of the subhalo mass function. 

Next to the results obtained using the concentration- 



FIG. 2. Boost using stripped-subhalo luminosities over that 
from field-lialo luminosities. Four fiducial models of the 
subhalo mass function are adopted, with minimum subhalo 
masses of m m i n = 10 -6 Mg and 10 4 Mg and slopes a = 1.9 
and 2. Starred symbols show results using the concentration- 
mass relation from Ref. [26J, assuming m m \n = 10 4 Mg and 
c v ir(^|m) oc z~ 0 ' 5 . 

mass relation of Ref. [|| (shown as solid and dashed 
curves in Fig. ED, we also show results for Milky-Way¬ 
sized halos when using the concentration-mass relation 
from Ref. [26j assuming m m i n = 1O 4 M 0 and c v i r (z|m) oc 
z~ 0 5 as starrled symbols. Both concentration models 
agree well for large mass halos, but differ significantly for 
smaller masses, closer to the resolution of the current- 
generation simulations, 10 5 M 0 . However, our results 
show that the boost ratio is insensitive to the initial 
choice of the concentration-mass relation. We also see 
that our semi-analytic model provides relatively smaller 
boost ratios compared with what is inferred from simu¬ 
lations directly [15j], as estimated above. It might be an 
indication that our approach provides a more conserva¬ 
tive boost relative to the dark-matter-only simulations, 
even though an increase in boost by upto a factor of ~4 
for the Milky-Way-sized halo is substantial. 

B. Boost 

Figure [3] shows the overall boost factor using the sub¬ 
halo mass functions that came out of our analysis (Ta¬ 
ble HD , as well as a few other phenomenological models of 
mass functions. For all cases we adopt the luminosities 
for stripped subhalos (solid lines), and compare to the 
luminosity using the ordinary field-halo approach (dot¬ 
ted). We caution that this boost can only be compared 
to other boost factors presented in the literature when 
taking the differences in the subhalo mass functions and 
concentration-mass relations properly into account. For 
example our boosts (solid lines in Fig. [3D are compara- 
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FIG. 3. Boost factors for halos of different mass using the 
concentration-mass relation from Ref. [§]. Solid curves in¬ 
clude the effect of tidal stripping, dotted curves assume field 
halo concentrations. The boost for three different subhalo 
mass functions are shown, using those from Table Q] (blue) 
and Ref. 0 (green and red). The expected boost for dwarf 
satellites of the Milky Way, adopting the mass functions from 
Table [I] is also shown (magenta). 


ble to those of Ref. EH where the tidal effect was not 
included. This is because our model is based on the 
concentration-mass relation from Ref. t 9|, which yields 
an even more modest boost and cancels the enhancement 
due to inclusion of the tidal effect. To be explicit, if we 
instead ran our analysis with the concentration-mass re¬ 
lation from Ref. 0 , we would have found a boost that is 
2-5 times larger than theirs. Similar arguments hold for 
different concentrations (including a simple power law). 


Estimate for dwarf spheroidal galaxies 

We estimate the expected boost for Milky-Way satel¬ 
lite galaxies. The density profile of the dwarf galaxies is 
taken to be that of a subhalo of given mass in a 10 12 M 0 
host. Therefore, the smooth component of the dwarf has 
a higher luminosity than that of similar-mass halos in 
the field. By de-projecting the surface brightness from 
substructures [35lT371 | , we estimate that about two thirds 
of the sub-subhalos lies outside of the tidal radius and is 
stripped away. This simple rescaling of the substructure 
mass function agrees with what is done by Refs. 0,0. 
However, this method likely yields an upper limit to the 
amount of sub-substructure, since, whereas sub-subhalos 
lose mass due to tidal effects, no additional sub-subhalos 
fall into the subhalo anymore [l5). The combined effect of 
the satellite being brighter than similar-mass field halos 
and the loss of sub-substructure makes the boost of satel¬ 
lite galaxies one order of magnitude smaller compared to 
their companions in the field. This supports the usual as¬ 


sumption that the boost due to sub-substructure is neg¬ 
ligible. Nevertheless, we show an estimate of how sub¬ 
substructure impacts our results in Fig. [3] For this esti¬ 
mate we assumed that two-thirds of the sub-substructure 
gets stripped away and that L ss h(m) = L s h(m). We ex¬ 
plicitly checked this at all host halo masses considered 
and the sub-substructure contribution is never more than 
- 10 %. 


VI. DISCUSSION 


We find that consistently modeling the subhalo lumi¬ 
nosity by taking into account tidal effects significantly 
enhances the global boost factor, compared to orthodox 
use of the concentration-mass relation. This result is in¬ 
dependent of uncertainties in the subhalo mass function 
or concentration-mass relation. 

Thus far, we applied a dark-matter only analysis, but 
state-of-the-art numerical simulations study the effects 
of baryons. Although they can change subhalo abun¬ 
dance and density profile, we do not expect them to have 
major impact on our results. First, the concentration- 
mass relations remain similar |39j, |40j . Second, low-mass 
(< 10 8 -10 9 Mq) halos, which give major contribution to 
the boost (Fig. [I]), are not expected to have a large bary- 
onic component in them. Nevertheless, we took a con¬ 
servative approach by estimating the boost ratio assum¬ 
ing that baryons would undo the effect of stripping com¬ 
pletely in subhalos < 10 8 Mq, and in the scenario where 
this has most impact (m m i n = 10 4 M 0 , and a = 1.9), the 
decrease is at most —30%. 

Encounters of subhalos with stars in the disk of the 
host will disrupt subhalos (e.g., Refs. [UHU). However, 
this happens only in a small volume close the halo center, 
and thus, will not affect the conclusions either. 

This study will have a broad impact on indirect dark 
matter searches in the extragalactic gamma-ray sky. 
Recent developments include the updated analysis of 
constraints on annihilation cross section from the diffuse 
gamma-ray background [43[, its anisotropies [371 44j, 
and cross correlations with dark matter tracers [45B48I ] . 
All these probes are subject to uncertainties in the 
halo substructure boost. Our conclusions are promising 
because having the boost factor larger by a factor of 
2-5 enhances the detectability (or improves the present 
upper limits) by the same factor. 
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APPENDIX 


First, we provide some details on the Vm ax -r max relation we find in our analysis and show that it is consistent with 
the one obtained with numerical simulations. Next, we have an extended discussion on the use of an Einasto profile 
rather than NFW. We repeat the order-of-magnitude estimate for the luminosity ratio in this context. Finally, we 
show how the boost depends on the minimum subhalo mass. 


A . 1 max r max relation 


In addition to the above analysis in terms of the evolution of the concentration parameters, we discuss our results 
in terms of the Wiax-fmax relation, whose evolution we modelled following Refs. [27[ as described above. 

For subhalos in the Milky-Way-sized host, we compare the I4iax“Imax resulting from our analysis to that of field 
halos with the concentration of Ref. Q. Like the simulation results [l5[T25j. we find that for the same Vm ax , subhalos 
have a smaller r max compared to field halos by a factor ^0.6. However, we find a softer slope (~1.13-1.15), which 
is a consequence of our choice for the concentration mass relation. Across seventeen orders of magnitude, we find 
(rmax,sh /rmax,fh ) ~ 0.6(Vmax,sh/Fmax,fh) and TTtsh ~ 6.2 X 10 7 (Vmax,sh/10 km S ) 3 AIq tO hold, leading to -b s h/Lfh ~ 
4. The I4iax-r max relation for our subhalos is plotted in Fig. IA-41 We also plot the relation for subhalos that is 
found in the Aquarius simulatio n [1511 . In addition, we show what is deduced for field halos when using the mass- 
concentration relation from Ref. [26l |. By running the analysis with this concentration-mass relation instead, we find 
a relation with a steeper slope that is consistent with the findings of the simulations [l]|. Concretely, we then find 
(Dnax,sh/Dnax,fh ) ~ 0.5(V^ na x,sh/^ / max,fh) 15 and TO s h « 4.9 x 10 7 (14iax,sh/10kms 1 ) 3,4 M 0 
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FIG. A-4. Knax-r m ax resulting from our analysis (dots) compared to that of field halos with the concentration of Ref. 
(green solid). The relation for field halos with the concentration of Ref. [26] (blue dotted) and the results from the Aquarius 
simulation[15|] (red dashed) are also shown for comparison. 


B. Einasto Profile 

Knax and r max are measurable quantities in the numerical simulations, and unlike the NFW scale radius and density, 
they are not profile dependent. However, in the above analysis we explicitly calculated Umax and r max starting from 
the assumption of an NFW profile for field and subhalos, and thereby we introduced a bias. Unfortunately, in our 
semi-analytic framework we are forced to resort to halo density profiles, only in simulations one is able to compare 
the U max -r max relation independently of the profile fl5| . 

The fact that our results resemble what is observed in simulations (Fig. IA-41) is encouraging. Nevertheless, we here 
also discuss what happens when applying an Einasto profile [23], |49| : 

/OEin(r) = p s ,Einexp 

Reference [5(j points out that, especially at large halo masses (> 10 14 Mq), the Einasto profile performs better than 
the NFW in fitting simulated halos. 

Below we will perform an order of magnitude estimate similar to that in main text, but now for an Einasto profile. 
In what remains we will closely follow the approach for calculating halo concentrations for Einasto profiles as laid out 
in Ref. (50j |. First, we define the concentration for the field halos with the Einasto profile as CEin = r v j r /r Si Ein- Starting 
from a concentration-mass relations for NFW halos, CEin can be obtained by requiring either that r maXj NFW = r ma x,Ein 
or Umax, nfw = Un aXi Ein- We will use the former. Since these are physical quantities, they should in principle be the 
same. However, we systematically deviate from the real U max and r max because we assume some density profile. As 
a result, by fixing r maXj NFW = 'Cmax.Ein we in general will obtain Vmax.NFW ^ Vmax.Ein- I n this case we obtain 

3.15 exp (—0.64a 1 / 3 ) 

CEin » cnfw- 2163 - ’ 8 

where we used r maXiE in ~ 3.15exp (—0.64a 1 / 3 ) r S) Ein (50j| . We obtain a from Eq. (23) in Ref. [5(j, which expresses a 
in terms of u, the peak height in the linear density fluctuation field. 

With this in place, we can now calculate the luminosity ratio as in the order-of-magnitude estimate in the main text. 
We again apply the concentration-mass relations from Ref. [26[ and use the phenomenological relations from Ref. [If| , 
which we repeat here: (r maXiSh /r maX)fh ) « 0.62(U max>sh /U maXifh ) 1 - 49 and m sh « 3.37 x 10 7 (U maX)Sh /10kms -1 ) 3 - 49 M©. 
We find CEin ~ 107 and Vmax.Ein ~ 1.3 km s -1 for rrifh = 10 5 . Finally, for the subhalos and field halos of equal mass 
with the Einasto profile, the relation, -Zf s h /Lfh ~ (Ps,sh/Ps,fh) (^s,sh/Cs,fh) = (bmax,sh/Unax,fh) (Cmax,fh/^max,sh); still 
holds. Note that the last equality is only exact if the field and subhalo have identical a’s. This yields L s h/Tfh « 4 for 
TO s h = 10 5 Mq and i s h/Tfh w 3 for m s h = 10 9 M 0 . 


^s.Ein 


- 1 


(7) 
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The above results are slightly lower than what we obtained when assuming an NFW profile. It can be understood 
by the fact that at higher masses Vmax.NFW < Vmax.Ein, when assuming r maXj NFW = Gnax.Ein- This makes field halos 
with an Einasto profile brighter than those with an NFW profile. However, most importantly, a change in profile does 
not alter our qualitative conclusions about the effect of tidal stripping on the subhalo luminosity. 


C. Boost dependence on mass of the smallest subhalos 

Finally, in Fig. lA-5l we show how the boost depends on the minimal subhalo mass in a MW-sized host. We show the 
boost for stripped subhalos (solid) and compare it to the field-halo approach (dotted) and again the concentration- 
mass relation is from Ref. _9|. This relation flattens at lower masses. We also show the power-law extrapolation from 
the Aquarius simulation [6|. This is equivalent to a power-law concentration-mass relation. 



lo s (w) 


FIG. A-5. The boost as a function of the minimum subhalo mass in a Milky-Way-sized host halo, for stripped subhalos (solid 
and dashed) and field-halo concentrations (dotted). The concentration-mass relation is from Ref. [9|, and flattens towards lower 
masses. We also show the power-law extrapolation of the Aquarius results [gj (red). 













